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The finite temperature behavior of small Silicon (Siio, Siis, and Si2o) and Tin (Snio and Sn2o) 
clusters is studied using isokinetic Born-Oppenheimer molecular dynamics. The lowest equilibrium 
structures of all the clusters are built upon a highly stable tricapped trigonal prism unit which is 
seen to play a crucial role in the finite temperature behavior of these clusters. Thermodynamics 
of small tin clusters (Snio and Sn2o) is revisited in light of the recent experiments on tin clusters 
of sizes 18-21 [G. A. Breaux et. al. Phys. Rev. B 71 073410 (2005)]. We have calculated heat 
capacities using multiple histogram technique for Siio, Snio and Siis clusters. Our calculated specific 
heat curves have a main peak around 2300 K and 2200 K for Siio and Snio clusters respectively. 
However, various other melting indicators such as root mean square bond length fiuctuations, mean 
square displacements show that diffusive motion of atoms within the cluster begins around 650 K. 
The finite temperature behavior of Siio and Snio is dominated by isomerization and it is rather 
difficult to discern the temperature range for transition region. On the other hand, Siis does show a 
liquid like behavior over a short temperature range followed by the fragmentation observed around 
1800 K. Finite temperature behavior of Si2o and Sn2o show that these clusters do not melt but 
fragment around 1200 K and 650 K respectively. 

PACS numbers; 61.46.+W, 36.40.-c, 36.40.Cg, 36.40.Ei 



I. INTRODUCTION 

The finite temperature behavior of small clusters is of 
great interest owing to the recent exciting calorimetric 
measurements carried out on clusters of sodium, tin, gal- 
lium and aluminumiiSi^*^. Calorimetric experiments on 
tin and gallium clusters with 10-50 atoms reveal many 
interesting features. For example, contrary to the stan- 
dard paradigm, it has been shown that clusters can have 
higher than bulk melting temperaturesSi^. Further, ad- 
dition of a single atom is seen to change the specific heat 
curve characteristics dramatically by changing a 'magic 
melter' into a so called 'non-melter— . These experiments 
have motivated many researchers to simulate the finite 
temperature behavior of small clustergSiSiSiifliiiiiSii^. 

In the present work, we study the finite tem- 
perature behavior of small silicon and tin clusters. 
The silicon clusters are of potential relevance to 
nanoelectronics industry and hence remain subject 
of many experimental as well as theoretical stud- 
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theoretical studies are devoted to the investigation of 
lowest equilibrium geometries of small silicon clusters, 
motivated by the experimental findings that Si clusters 
undergo structural transition (prolate to spherical) in the 
range 22-34Mi2&. To date, global minima of silicon clus- 
ters with sizes n<24 have been well established from var- 
ious recent unbiased searches such as genetic algorithm, 
single-parent evolution algorithm, 'big bang' optimiza- 



tion et(}iSiS2*2&. 

The surface induced dissociation studies revealed that 
the fragmentation behavior is common within the clus- 
ters of semiconducting group 14 elementa^SiS. Studies 
on silicon clusters of various sizes report existence of 
fragmentation pathway as compared to evaporation pro- 
(,gggi5j3i predict that clusters with sizes 6-10 are the 
most abundant fragments and are referred to as magic 
fragmentsSii^. Although, it is speculated that silicon 
clusters up to sizes 70 will undergo fragmentation rather 
than evaporationi^, it is not yet clear if the clusters un- 
dergo the traditional solid-like to liquid-like transition 
before fragmenting or they fragment without melting as 
reported recently in case of tin clusters'^^. Hence, we 
report detailed thermodynamic simulations on Si„ clus- 
ters with various sizes; n — 10, 15, 20. These simulations 
have been carried out using Density Functional Theory 
(DFT) within the Generalized Gradient Approximation 
(GGA). The total simulation time for each cluster is at 
least 2 ns. Our detailed analysis of finite temperature 
behavior shows that the finite temperature behavior of 
Siio is dominated with isomerization and it is rather dif- 
ficult to identify the region corresponding to solid-like to 
liquid-like transition. However, the cluster clearly frag- 
ments around 2800 K. Sii5, on the other hand, does ex- 
hibit a liquid like phase over a short temperature range 
before fragmenting around 1800 K. Si2o does not melt 
but fragments around 1200 K. 

Both silicon and tin are group IV elements and it is 
interesting to note some peculiarities of this group. Car- 
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bon, the first element is a non-metal with a high energy 
gap between valence band and conduction band (around 
5.5 eV) whereas silicon and germanium are semiconduc- 
tors at room temperature. Tin and lead are metals at 
room temperature of which tin undergoes a structural 
phase transition and transforms into a semiconductor be- 
low 286 K. Contrary to the bulk, tin clusters (n<25) re- 
semble greatly to that of small silicon and germanium 
clustersSi^. Interestingly, Si, Ge and Sn clusters with 
sizes > 9 are built on stable TTP unit which, as we shall 
see, plays a crucial role in finite temperature behavior of 
these clusters. However, in contrast to silicon clusters, 
tin clusters undergo a size rearrangement from prolate 
(stacked TTP units) to a more compact spherical shape 
in a considerably broad range (from 35-65 atom cluster). 
The ionic mobility experiments exploited the fact that 
small tin clusters have prolate ground states. Experi- 
ments on small tin clusters (10-30 atom clusters) indi- 
cated that these clusters do not melt at least 50 K above 
the bulk melting temperaturei^ This conclusion is based 
on the argument that presence of a liquid phase is identi- 
fied with the change in the shape from prolate to spheri- 
cal and hence enhanced ionic mobilities in the liquid like 
regior^. Density functional simulations supported these 
findings showing that, indeed, melting temperatures of 
small tin clusters are at least 1000 K higher than that of 
Trrafbttifc ]^'^'^'"*"^ - Our previous calculations, within LDA 
showed that Snio has substantially higher melting tem- 
perature, 2300 K whereas Sn2o is in a liquid-like phase 
after 1200K§^. However, the recent experiments on Sn^g, 
Snfg, SujQ and Suj^ demonstrated that these clusters 
do not melt, but sublimate^ around 650 K. In light of 
these findings, we simulate the thermodynamics of tin 
clusters, Snio and Sn2o, using GGA functionals. In the 
present work, simulations are extended over much longer 
time scale, at least 90 ps per temperature as compared 
to our earlier simulations of 40 ps per temperatureSi^. 
Our studies confirm the fragmentation observed experi- 
mentally and bring out similarities between small tin and 
silicon clusters in terms of structure and dynamics. 

In what follows, we present computational details in 
Sec. m Results of both tin and silicon clusters are given 
in Sec, mil We discuss various issues concerning the struc- 
ture and dynamics of all the clusters in the same section. 
We conclude our paper in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

All the thermodynamic simulations are performed us- 
ing Born-Oppenheimer molecular dynamics based on 
Kohn-Sham formulation of DFTi^l. We have used Van- 
derbilt's ultrasoft pseudopotentiala^ within GGA, as im- 
plemented in VASP package^ for all the clusters. Ther- 
modynamic behavior of tin clusters is also studied using 
LDA as implemented in VASP. For all the calculations, 
we use only 3s^ and 3p^-electrons as valence in case of Si. 
In case of Sn, we use 5s^ and 5p^-electrons as valence. 



taking d-electrons as a part of the ionic core. An energy 
cutoff of 13.84 Ry and 9.77 Ry is used for the plane-wave 
expansion of Si and Sn wave functions respectively, with 
a convergence in the total energy of the order of 10~4 cV. 
Cubic supercells of lengths 15, 18 and 22 A are used for 
Xio, Xi5 and X20, where X — Si,Sn, respectively. 

The ground state and other equilibrium structures of 
all the clusters are found by optimizing several struc- 
tures chosen periodically from a high temperature sim- 
ulation. We analyze the bonding characteristics of the 
ground state and few excited state isomers of the clus- 
ters using Electron Localization Function (ELF)i2i. This 
function is normalized between zero and unity; a value of 
1 represents a perfect localization of the valence charge 
while the value for the uniform electron gas is 0.5. The 
locations of maxima of this function are called attractors, 
since other points in space can be connected to them by 
maximum gradient paths. The set of all such points in 
space that are attracted by a maximum is defined to be 
the basin of that attractor. Basin formations are usually 
observed as the value of the ELF is lowered from its max- 
imum, at which there are as many basins as the number 
of atoms in the system. Typically, existence of an iso- 
surface or a basin along the bonding region between two 
atoms at a high ELF value, say > 0.7, signifies a localized 
bond in that region. 

For examining the finite temperature behavior, ionic 
phase space of all the clusters is sampled by isokinetic 
MD where kinetic energy is held constant using velocity 
scaling. In case of Xio {X = Si, Sn) clusters we split the 
total temperature range from 100-3000 K into 22 different 
temperatures. In case of Sii5 and Si2o clusters we split 
the temperature range from 100-1800 K into 14 different 
temperatures. Each system is heated to desired temper- 
ature from previous temperature slowly within a time 
scale of 15 ps. Molecular dynamics is then simulated at 
each temperature for 30 ps after which it is subsequently 
heated to the next temperature. The MD for each tem- 
perature is then continued for an additional time scale of 
90 ps so as to have sufficiently large statistics. In case 
of Sn2o , the molecular dynamics simulations were carried 
out within the LDA as well as GGA approximations. In 
case of GGA, we have studied the finite temperature be- 
havior of the cluster around 200 K, 350 K, 500 K and 
650 K. In case of LDA calculations the temperature range 
is split so to have the finite temperature behavior of the 
cluster around 100 K, 250 K, 500 K, 650 K, 800 K and 
1000 K. 

Following the MD simulation, classical ionic density 
of states of system is extracted using multiple his- 
togram technique and the details are found in litera- 
ture^^. Various other thermodynamic indicators such as 
Mean Square Displacements (MSD) of atoms and Root- 
Mean Square bond-length fluctuations (RMS-BLF, <5rms) 
are computed. More technical details concerning the 
extraction of thermodynamics averages, indicators and 
computation of specific heat curve can be found in our 
earlier paper—. 
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FIG. 1: The ground state, and some representative excited 
states of Siio, Siis and Si20- The energy difference Ai5 is 
given in eV with respect to the ground state. 



III. RESULTS AND DISCUSSION 
A. Thermodynamics of Silicon clusters 

First, we discuss the finite temperature behavior of 
Siio- The ground state and various excited state geome- 
tries of Siio are shown in FigQ]- (a) to (d). The ground 
state of Siio is a tetra capped trigonal prism (similar to 
that of Snio^) and is consistent with the one reported 
earlier—. It may be noted with some interest that the 
ground state and the first excited state structure (higher 
by 0.2 eV) contain one TTP unit whereas all other ex- 
cited configurations do not contain a TTP unit. This 
explains the significant energy difference (of « 1 eV) be- 
tween the first two equilibrium structures and the re- 
maining excited state structures. We also note that in 
the ground state geometry, all the atoms except one cap- 
ping the triangular face of trigonal prism have a mini- 
mum of four-fold coordination. The nature of bonding is 
analyzed using ELF which shows the expected covalent 
bonding among all the silicon atoms. 

Next, we discuss the calculated ionic specific heat curve 
for Siio shown in Fig.|21 The salient features of the curve 
are presence of a broad shoulder from 350 K to 1400 K 
and a main peak at 2300 K. An analysis of ionic mo- 
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FIG. 2: Normalized canonical specific heat of Siio. 



tion at various temperatures, reveals several additional 
aspects. The shoulder at 350 K arises from a structural 
transition of cluster from the ground state to the first 
excited state and back to the ground state. This process 
is due to the attempt of the tri-coordinated atom (the 
atom capping the triangular face of the trigonal prism) 
to acquire the forth neighbor. The isomerization seen 
around 350 K is quite similar to that of Snio'* and occurs 
in such a way that at least three atoms (cap on triangular 
face, rectangular face and one from the trigonal prism) 
exchange their positions without changing the geometry 
considerably. This transition occurs with increased fre- 
quency until 1400 K and leads to the diffusion of atoms 
through out the cluster without considerable change in 
the geometry. It is only at 1600 K and beyond that struc- 
tures corresponding to the other excited states (which 
are higher in energy by almost 1 eV) are observed. Thus, 
the rise in the specific heat curve observed after 1400 K 
is associated with the destruction of the TTP unit and 
occurrence of other excited states. At much higher tem- 
peratures, around 2300 K and above, we observe sev- 
eral high lying configurations which can be thought of as 
two weakly bonded clusters of smaller sizes, with ener- 
gies above 2.5 eV with respect to the ground state. The 
cluster eventually fragments around 2800 K into Sir and 
Sia as shown in Fig. 01 

Thus, though the main peak in the specific heat curve 
of Siio is seen around 2300 K, it cannot be attributed to 
the solid- like to liquid-like transition. In fact, we believe 
that this peak is due to the fragmentation that occurs 
around 2800 K. This is also evident from the analysis 
of the root mean square bond-length fluctuations (<5rrns) 
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FIG. 3: Siio seen as fragmenting into Sir and Sis at 2800K 



Siio 






(a) 



0.5 


Sn,o 




0.4 






CO 






0.2 
0.1 




(b) 









500 1000 1500 2000 2500 3000 

Temperature (K) 



500 1000 ISOO 2000 2S00 3000 

Temperature (K) 



sii5 ; 


0.75 
0.<i 




« 0.45 




{).3 


/ (c) 


0.15 




400 800 1200 1600 20tK) 

Temperature (K) 



300 600 9(H) 1200 1500 

Temperature (K) 



FIG. 4: Root mean square bond length fluctuations (5ri 
for Siio, Snio, Siis and Si2o. 



shown in Fig. 21"(a). It is clearly seen that 6^^^ is char- 
acterized by two sharp increments, first one from 175- 
650 K and another at 2800 K. According to the bulk 
Lindemann criterion a value of (5rms greater than 0.1 is 
indicative of transition from solid phase to liquid phase. 
In bulk, the solid-liquid transition is sharp whereas in 
clusters, the transition is gradual and is observed that, 
in a liquid-like phase, the value of (5rms saturates around 
0.2 - 0.3. At 650 K J^ms reaches a value of 0.2. This 
is due to the fact that GS-First Excited state-GS transi- 
tion results in diffusion of atoms through out the cluster 
thereby saturating the value of (5rms. The value of Jrms 
is almost constant for next 1000 K (i.e. up to 1600 K). 
Above 1600 K, coincident with the destruction of TTP 
unit, the Jrms experiences a gradual and slow rise until 
2600 K. This gradual rise ends with a second sharp jump 
at 2800 K due to fragmentation of the cluster. 
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FIG. 5; Atomic distribution from the centre of mass for Siio 



Next we analyze the average distribution of atoms from 
Center of Mass (COM) of the cluster as a function of tem- 
perature (shown in Fig. . We have plotted the proba- 
bility of finding an atom as a function of distance from 
the COM. At low temperatures, a perfect shell structure 
is retained, as observed for 100 K (shown in Fig. El-a). 
It can be easily identified that the first two peaks corre- 
spond to the two faces of the trigonal prism, the third 
peak corresponds to the caps on the rectangular surfaces 
and the last one is due to the cap on the triangular face. 
With increasing temperature, the peaks begin to merge 
in the bottom, (as seen around and above 350 K) indicat- 
ing an exchange of atoms between various shells. Around 
650 K, four peaks merge into two and the distribution of 
atoms remains identical till 1400 K. Above 1600 K, the 
COM plot indicates absence of any sort of average order 
in the arrangement of the atoms. The peak gets more 
broader at higher temperatures and around 2300 K (co- 
incident with the peak in the specific heat), the atomic 
distribution from COM is seen as a broad curve. 

In short, the main peak in the specific heat curve oc- 
curs around 2300 K which is due to fragmentation. The 
Lindemann criteria, atomic distribution of atoms around 
COM and other qualitative parameters indicate diffusive 
motion of atoms within the cluster around 650 K. Al- 
though, for a very broad range of temperatures, from 
650 K to 1400 K, the cluster undergoes characteristic iso- 
merization discussed earlier, it retains the shape of the 
cluster, but allows diffusion of the atoms through out the 
cluster. Thus, finite temperature behavior of very small 
clusters like Siio, is dominated by isomerization. The 
solid-like to liquid-like transition is continuous and it is 
rather difficult to discern the transition temperatures for 
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the same. This is also confirmed by our on going studies 
on small sodium and lithium clusters^. 

Turning to the finite temperature study of Siis, we 
first discuss the equilibrium structures of Siis shown in 
Fig. ^-(e to h). The ground state geometry of Siis can 
be described as a tricapped trigonal prism fused with a 
tricapped trigonal anti prism (Cs^ symmetry Fig.^(e)). 
This is consistent with the lowest energy configuration 
reported by CCSD(T)2? and DFT simulations^^ . The 
first excited state consists of two TTP units sharing a 
common triangular face as shown in Fig. ^-(f). This 
structure can be viewed as a serious distortion of the 
lower six atom ring present below the TTP unit in the 
ground state. The first few excited states are built on 
TTP unit of Sig as shown in Fig.^(f,g), whereas isomers 
observed at higher temperatures can be thought of as a 
combination of two small silicon clusters as shown in Fig. 
^-(h) and suggest probable paths for fragmentation. 

The ionic specific heat curve computed for Siis (see 
Fig. has two small peaks around 400 K and 800 K and 
a main peak around 1400 K which is associated with frag- 
mentation of the cluster around 1800 K. Ionic motions of 
Siis cluster at various temperatures provide a plausible 
explanation for different features observed in the specific 
heat curve. We would also like to note that Siis has many 
low lying isomers (we have at least 30 different equilib- 
rium structures out of which few relevant structures are 
shown in the figure) unlike Siio- The isomerization be- 
gins around 400 K when we observe the first excited state 
within our calculation (Fig. ^-(f)). This involves rear- 
rangement of lower six atom ring. With further rise in 
temperature, the cluster remains in the first excited state 



for a reasonable amount of time after which the atoms 
in upper TTP unit reorient into a six atom ring leading 
to an inverted ground state. Thus, the first excited state 
appears to be an intermediate for interchanging the posi- 
tions of TTP unit and the hexagonal ring in the ground 
state. This isomerization process does not involve diffu- 
sion of atoms within the cluster and continues to occur 
with increased frequency until 700 K. Further rise in tem- 
perature leads to several other isomerization routes along 
with the one described above which results in to diffusion 
of atoms through out the cluster. This corresponds to the 
second shoulder seen in the specific heat curve around 
800 K (see Fig. Around 1200 K isomers correspond- 
ing to two small silicon clusters (Fig.^(h)) are observed 
in the ionic motion which provide a fragmentation path 
at higher temperatures. 

The (Sj-ms (shown in Fig.2|-(c)) and the distribution of 
atoms from COM of Siis (shown in Fig.O support these 
observations. Linear rise in S^ms at lower temperatures is 
indicative of pure vibrational motion. Beginning of iso- 
merization is shown up as a rise in drms around 400 K 
which continues till 900 K. The value of S^ms saturates 
after 900 K showing a liquid-like state over a short tem- 
perature range. Around this temperature, the diffusive 
motion of atoms between the upper and lower units of 
Siis is clearly seen in the distribution of atoms from the 
COM plot. The value of i5rms increases abruptly around 
1800 K due to the fragmentation of Siis into Sig and Sig 
clusters. 

It is interesting to compare the finite temperature be- 
havior of Siis and Siio. Although, both clusters frag- 
ment, Siis (1800 K) dissociates at considerably lower 
temperature than Siio (2600 K)^^. In case of Siio, all 
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isomers except one, are at much higher in energy. As a 
resuh, the cluster undergoes a pecuUar isomerization (de- 
scribed earher) over a wide range of temperatures (350 K 
to 1400 K) leading to a broad shoulder in the specific 
heat curve and a flat region in the Jrms plot. On the 
other hand, in case of Sii5, the isomerization region is 
spread over comparatively shorter range of temperatures 
(400 K to 900 K). 

As we have already noted, for the clusters of group 
IV elements, TTP forms a basic block on which larger 
clusters are built upon. As will be discussed, the ground 
state of Si2o is a prolate structure whereas the TTP unit 
is evident in the excited states. Though the fragmen- 
tation of Si2o is expected, it will be interesting to see 
if the fragmentation is preceded by a liquid like state. 
In what follows, we present the thermodynamics of Si2o 
cluster. In Fig. ^"(i) to (1), we show the ground state 
as well as some of the excited state geometries of Si2o- 
The lowest-energy isomer of Si2o is composed of three 
units: a so called magic number cluster Sig, a hexagonal 
chair unit in the middle, and a low energy isomer of Sis 
unit. The first excited state consist of two distinct TTP 
units whereas the second one is non-prolate, spherical 
structure. The forth equilibrium structure that we have 
shown in the figure, (see Fig. n]-(l)) has two Siio units 
connected via very weak bonds (bond lengths « 2.7 A 
as compared to the normal Si-Si covalent bondlength of 
2.34 A). As we shall see, the cluster fragments via this 
isomer. All the equilibrium structures shown here are 
consistent with earlier reportsSiS^. 

The nature of bonding is examined using ELF. As ex- 
pected, ELF for the ground state shows a high covalent 
bonding with all the basins completely connected by an 
isovalue of 0.75. It is also observed from ELF that the two 
TTP units in Fig.^(l) are weakly bonded to each other. 
Turning to the finite temperature behavior of Si2o, we 
note that till 800 K, atoms in the cluster vibrate about 
their mean positions and no isomerization is observed. 
Around 1000 K, the cluster transforms from the ground 
state to the excited state having two distinct Siio units 
(Fig-CKl)) and eventually fragments into two Siio units 
around 1200 K. This is refiected in the S^ms shown in 
Fig. 0]-(d) where the fragmentation is evident from the 
sharp rise around 1200 K. In other words, finite temper- 
ature behavior of Si2o does not show solid-liquid transi- 
tion prior to fragmentation. This may be attributed to 
the low dissociation energy barrier in case of Si2o cluster 
(of the order of 1.2 eV) in comparison to that of Siio and 
Sii5 (which have a dissociation energy barrier of 4.2 eV 
and 2.2 eV respectively)^^. To understand the infiuence 
of starting geometry on the fragmentation behavior of 
the cluster, we have repeated the thermodynamic sim- 
ulations on Si2o cluster with Fig. ^-(j) and Fig. n]~(k) 
as the starting configurations. We observe the cluster 
to fragment around 1200 K in both these calculations 
as well indicating that the fragmentation behavior is not 
sensitive to the starting configuration of the cluster. 

We end our discussion on silicon clusters with an inter- 





FIG. 8; Prolate and spherical configurations of 8120 around 
1800 K 



esting observation. The observed fragmentation can be 
bypassed with higher rate of heating. In all these simu- 
lations, the cluster is heated to a particular temperature 
from the previous temperature at a rate of 1 K in 30 fs 
and then maintained at that temperature for next 30 ps. 
This is used as a starting point for heating the cluster 
to the next temperature. We observed the fragmenta- 
tion after maintaining the cluster at 1200 K for at least 
60 ps. On the other hand when the cluster is heated 
from 1200 K to 1600 K, without sufficiently thermaliz- 
ing it at 1200 K, it shows liquid-like behavior at 1600 K. 
This liquid like state is also observed for 1800 K^^. It 
is seen that while heating the cluster without allowing 
enough time for thermalization, it attains enough kinetic 
energy to push the two Siio units closer to each other 
(SiTTPi-SiTTP2 bond distance of 2.2 A). Consequently, 
the two TTP units begin to interact leading to high ly- 
ing configurations (containing no TTP units) as shown 
in Fig. IHl The ionic motion of this super heated cluster 
shows that the cluster does not fragment at the end of 
100 ps. On the contrary, a complete diffusive motion of 
atoms within the cluster is observed. The temperature of 
the cluster was further raised to 1800 K. The cluster was 
then maintained at that temperature for an additional 
time scale of 100 ps. An analysis of ionic motion at this 
temperature shows the cluster to be in a liquid-like state 
there by revealing a possibility to bypass the fragmenta- 
tion process. An examination of ionic motion at 1600 K 
and 1800 K reveals many interesting features. Around 
1600 K, the cluster oscillates between prolate configura- 
tions and non-prolate (near spherical) configurations sev- 
eral times. As the temperature rises, the cluster spends 
more time within the near spherical configurations. This 
shape change can be experimentally detected using ion 
mobility experiments (diffusion coefficient)^. Hence, is 
it interesting to analyze these shape changes in a more 
detailed fashion. In Fig.O we have plotted the deforma- 
tion coefficient (cpro) of the cluster as a function of time, 
at 1600 K and 1800 K respectively. For a given ionic 
configuration, Sp^o is defined as 
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FIG. 9: Deformation coefficient, epro, as a function of simu- 
lation time for Si2o cluster 



where, Qi,Q2 and Q3 are the eigenvalues, in descending 
order of the quadrapole tensor 

Qij = ^ RiiRij (2) 
/ 

I runs over the number of ions and R/i is the ith coor- 
dinate (i and j run from 1 to 3) of the ion 'I' relative to 
the center of mass of the cluster. An ionic configuration 
with spherical shape has a value of tpro — 1- Devia- 
tion of the value of from 1 is a measure of prolaticity 
of the cluster. As is seen from the figure, the cluster 
oscillates between prolate and nearly spherical configu- 
rations quite frequently around 1600 K. At 1800 K, the 
oscillations reduce and the structure is seen to visit the 
spherical configurations more often. The value of Jims 
at 1600 K and 1800 K is 0.35 and 0.37 respectively and 
satisfies the Lindemann criterion. 



B. Thermodynamics of Tin clusters 

Various interesting observations seen in case of the sil- 
icon clusters and recent experiments^ on tin clusters 
(Snig, Snig, Sn2o and Sn2i) have motivated us to revisit 
the thermodynamics of tin clusters (Snio and Sn2o). In 
the present work, the finite temperature study is carried 
out within the GGA as well as LDA approximations us- 
ing the vasp package^. The clusters are heated to each 
temperature at a much slower rate (1 K in 0.036 ps) than 
in our earlier work. The cluster was maintained at each 
temperature for a period of aleast 90 ps (as contrast to 



40 ps per temperature in our earlier work) . The heat ca- 
pacity curves are computed from the above calculations 
(GGA and LDA) and compared with our earlier results. 

The finite temperature behavior of Snio is seen to be 
considerably similar to that of SiiQ. This is in fact ex- 
pected result as the ground state and even the excited 
state configurations of Siio and Snio are identical. It is 
gratifying to note that the specific heat curves computed 
from GGA as well as LDA calculations are nearly simi- 
lar to that of our earlier reported one — (with a slightly 
shifted main peak around 2200 K as compared to our ear- 
lier reported one at 2300 K). Various inferences concern- 
ing bonding and ionic motion are in perfect agreement 
with our earlier report. The significant difference to be 
noted is the fragmentation of Snio observed at 2600 K. 
This process is observed only after a time scale of 60 ps. 
We believe this feature to be missing in our earlier re- 
ports due to the smaller simulation time scales ( 50 ps 
per temperature). However, it may be recalled, that ionic 
motion of cluster at 2600 K in our earlier study revealed 
few configurations corresponding to weakly interacting 
Sn^; and Snio_a; units. Thus, we conclude that the peak 
in the specific heat curve is due to the fragmentation of 
Snio cluster seen at 2600 K and not due to the traditional 
solid-liquid transition as reported. 

We have also computed the ionic specific heat of Sn2o 
cluster within GGA and LDA approximations. Analy- 
sis of the ionic motion of Sn2o cluster (GGA) at vari- 
ous temperatures shows that the cluster vibrates around 
its ground state until 600 K. Around 650 K, the cluster 
fragments into two TTP (Snio) units. This is in per- 
fect agreement with the recently reported experimental 
observations ~. Finite temperature study of Sn2o clus- 
ter within LDA approximation shows the cluster to frag- 
ment into two Suio units around 1000 K. The failure to 
detect the fragmentation behavior in our earlier reports 
is attributed to the following reasons: (a) faster heat- 
ing rate (1 K in 0.01 ps) with smaller simulation time 
scales in our earlier simulations (40 ps per temperature), 
(b) Use of LDA functional in our earlier reports. As 
already mentioned ground state Sn2G^ cluster has inho- 
mogeneous bond strength distribution. While the bonds 
within intra-TTP atoms are strongly covalent in nature, 
intcr-TTP atoms are connected through very weak cova- 
lent bonds. The LDA approximation was unable to dis- 
tinguish these differences leading to an over-estimation 
of binding energies between the two TTP units. Hence, 
the Sn2o cluster did not fragment in the small simulation 
times scales reported earlier. Snio cluster in contrast, 
has all its atoms connected with similar bond strengths, 
resulting in a more accurate finite temperature behavior 
within the LDA approximation. 



IV. CONCLUSIONS 

We have presented thermodynamics of small silicon 
and tin clusters (Siio, Sii5, Si2o, Snio and Sn2o) within 
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LDA and GGA approximations of DFT. Finite temper- 
ature behavior of Siio and Snio is dominated by iso- 
merization process and fragmentation is observed around 
2800 K and 2600 K respectively. The peak in the spe- 
cific heat curves around 2300 K and 2200 K for Siio and 
Snio correspond to the observed fragmentation rather 
than sohd-Hke to hquid-Hke transition. The similari- 
ties observed in the finite temperature behavior of these 
two clusters suggest the influence of ground state struc- 
tures on the finite temperature behavior of a cluster. 
The specific heat curve of Siis has a main peak around 
1400 K which is due to the fragmentation observed 
around 1600 K. However, the cluster exhibits a liquid 
like behavior over a short temperature range (900 K to 



1400 K) before fragmenting. Si2o and Sn2o are found 
to dissociate without melting around 1200 K and 650 K 
respectively. Our simulations on Sn2o agree with the re- 
cently reported experimental results. 
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